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Abstract 

Supernova remnants are usually analysed in the light 
of hydrodynamical models of the interaction of supernova 
ejecta with either a constant density ambient medium or 
a circumstellar medium produced by a constant presuper- 
nova wind. However, the ejection of energetic wind dur- 
ing the presupernova phase changes the ambient medium 
structure and, consequently, the early supernova remnant 
evolution. We have analysed the evolution of young rem- 
nants of type la supernovae, focusing on the imprint of 
the presupernova wind history on the supernova remnant 
structure and on the influence of the explosion mechanism. 
We have found that the remnant evolution is most sensi- 
tive to the explosion mechanism at ages not larger than a 
few hundred years, while the presupernova history shows 
its influence at later epochs, before the Sedov phase sets 



that allow more realistic and more detailed supernova 
(SN) models. Here, we have focused on the last kind of 
studies. 

Presupernova evolution is often neglected in the hy- 
drodynamic simulations of type la SNR, because the mass 
loss episodes of the progenitors are not believed to have a 
deep influence on their circumstellar medium (CSM) . The 
shaping of the CSM through stellar winds and the for- 
mation of wind bubbles have been studied extensively for 
massive early type stars (i.e., type II SN progenitors), al- 
beit only for constant mass loss episodes or combinations 
of them ( |Weaver et al. 1977| , |Chevalicr fe Liang 1989) ). Up 
to now, the usual assumptions for the CSM surrounding 
type la SN progenitors have been either a constant density 
interstellar medium (ISM) or a power law CSM merging 
smoothly with the ISM. However, recent models of type la 
presupernova evolution predict substantial non-constant 



mass loss episodes during the binary evolution (Hachisu 



O ; Key words- rirrumstellar matter - stars- winds, outflows ^ato fe Nomoto 1996|]Langer et al. 2000|). Evidence for this 



supernova remnants - supernovae: general 



1. Introduction 

Supernovae have become one key piece in the knowledge 
building of modern astronomy. They stand at the corner- 
stone of many fields, including cosmology, galactic and 
stellar evolution, cosmic rays generation, and many oth- 
ers. Type la supernovae (SNIa) are currently used as prac- 
tical tools to measure cosmological distances and obtain 
insight about what kind of Universe we are in. However, 
the theoretical knowledge of SNIa is still plagued by un- 
certainties about the nature of their progenitor systems 
and the nature of the explosion mechanism. 

Modern X-ray observatories, like Chandra and XMM- 
Newton give the opportunity to collect a huge amount of 
high quality, high resolution, data about supernova rem- 
nants (SNRs). Whereas the atomic codes that allow to 
compute and fit the X-ray spectra have experienced a 
considerable improvement in the last years, this is not 
the case of our knowledge of the hydrodynamical evolu- 
tion of SNRs. There are two desirable directions of ad- 
vance for our theoretical understanding of SNR evolution: 
one is multi-dimensional hydrodynamical simulations, and 
the other is one-dimensional hydrodynamical calculations, 



kind of behavior has been found in X-rays even for type II 



events (Emmler, Aschenbach & Wang 2001). Light echoes 
have been detected around two ovcrluminous type la su- 
pcrnovae: SN 1991T ([Schmidt et al. 1994Q and SN 1998bu 
(Cappellaro et al. 2001), which have been attributed to the 
pressence of circumstellar dust. Thus, it is now necessary 
to evaluate the influence of such non-constant presuper- 
nova winds on the properties of the SNR. 

The explosion mechanism of SNIa is still a matter 
of strong debate. Since the evolution of young SNRs is 
dominated by the SN ejecta profile, different explosion 
mechanisms (implying different ejecta profiles) can lead 
to imp ortant variations in the remnant behavior at early 
stages. Dwarkadas fc Chevalier (1998)| found that most 
ejecta profiles could be adjusted by an exponential pro- 
file, but also that this approach fails to reproduce many 
features of SNR evolution. Their results call for a more 
thorough study of the influence of the explosion mecha- 
nism on the properties of SNRs. 

2. Models and results 

We have followed with an hydrodynamical code the inter- 
action of the presupernova wind with the surrounding ISM 
up to the time of SN explosion . The code was the same 
as in Badenes fc Bravo (2001) , modified to take into ac- 



count the radiative losses, which are dynamically meaning- 
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Table 1. Charateristics of wind models. 







M w 


K 51 


^■cnd 


tsN 


R s 


^s, Weaver 


Wind model 


(km s- 1 ) 


(Mq) 


(xlO 51 ergs) 


(Myr) 


(Myr) 


(xlO 19 cm) 


(xlO 19 cm) 


A 


200 


0.2 


8 x 10~ 5 


0.2 


0.7 


3.2 


2.6 


B 


20 


0.2 


8 x 10" 7 


0.2 


0.7 


1.4 


1.0 


C 


200 


0.6 


2.4 x 10" 4 


1.5 


1.5 


5.0 


4.5 


D 


20 


0.6 


2.4 x 10" 6 


1.5 


1.5 


2.0 


1.8 



NOTE.- v„ is the wind velocity, M w is the total mass ejected in the wind, and K&i is the total kinetic energy of the wind; 
tend is the duration of the wind phase, and tsN is the time of the SN explosion; R s and -Bs.Wcavcr are the radius of the outer 
shell at the time of SN explosion, according to the present calculations and to Weaver et al. (1977), respectively. 
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Figure 1. Time evolution of wind mass-loss rate for mod- 
els A and B (dashed line) and models C and D (solid line). 
The time at which the SN explodes is identified by a star. 
Note that in models A and B the wind ends 0. 5 Myr before 
SN explosion. 



four cases corresponding to different wind models, whose 
characteristics can be found in Table [l] and Figure @ (the 
ISM was always assumed to be initially homogeneous). 
Figure § displays the CSM profile for each model at the 
time of the SN explosion. 

The results of the hydrodynamic computation of wind- 
ISM interaction have been used as the start configuration 
for the hydrodynamic calculation of SNR evolution. In this 
study of the influence of different wind models we have 
used an analytic ejecta profile given by an exponential 
law. The results are shown in Figures |^ and |], together 
with the evolution of the same ejecta profile interacting 
with a constant density ISM. We use the shock expansion 
parameters (77 = dlnr/dlni) of the forward and reverse 
shocks as a figure of merit for model comparison. It can 
be seen that the wind ejected by the progenitor system of 
a SNIa can have a dramatic influence on the dynamical 
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Figure 2. Circumstellar medium density profile for each 
wind model at the time of SN explosion. Lines correspond 
to wind models A (solid line), B (dotted line), C (dashed 
line), and D (dash-dotted line). 



evolution of its SNR. This evolution cannot be accurately 
described with the usual assumptions of a power law CSM 
merging smoothly with the ISM or direct interaction be- 



ful during the wmd-lSM interaction. We have computed 2001) 



tween the ejecta and the ISM (see also Badenes & Bravo 



We have also made a preliminary study of the influence 
of the explosion model on the SNR properties. We have 
chosen a few representative models of SNIa explosion (see 
Table |J): a delayed detonation (DDT), a pulsating de- 
layed detonation (PDD), a deflagration (DEF) and a sub- 
Chandrasekhar explosion (SCH) . In Figure ||, we compare 
them with the most usual analytical profiles: an exponen- 
tial (EXP) and a power law with index 7 (PWL), both 
with ejected mass 1.4 Mq and kinetic energy 10 51 erg. 

The interaction of the different explosion models with 
a homogeneous ISM has been followed with the same hy- 
drodynamical code. The evolution of the forward shock 
expansion parameter is given in Figures I and (the evo- 
lution of the reverse shock expansion parameter shows 
only minor deviations between the different models and 
is not shown). It can be seen that the dynamical behav- 
ior of SNRs resulting from realistic SNIa explosion models 
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Table 2. Characteristics of the SNIa explosion models. 



Explosion 


M ej 


K 51 


t' 




Mca 


M s 


M Si 






model 


(Mq) 


(xlO 51 ergs) 


(10 9 s) 


(Mq) 


(Mq) 


(Mq) 


(Mq) 


(Mq) 




DEF 


1.38 


0.79 


11.4 


0.68 


0.006 


0.022 


0.031 


4.7 x 10" 


-6 


DDT 


1.37 


2.00 


7.1 


1.04 


0.023 


0.069 


0.084 


5.6 x 10" 


-5 


PDD 


1.36 


1.52 


8.2 


0.83 


0.020 


0.070 


0.090 


1.1 x 10" 


-5 


SCH 


0.97 


0.98 


7.7 


0.49 


0.027 


0.102 


0.141 


1.2 x 10" 


-3 



NOTE.- M j and K51 are.respectively, the ejected mass and the total kinetic energy r eleased, and t' is th e normalization 
timescale (see Figs. ^|and|^). The rest of columns give the yields of selected species. See Bravo et aL (1996) for more details 



on the models. 




Figure 3. Expansion parameter of the forward shock as 
a function of time for the interaction of an exponential 
ejecta profile with the CSM as modifed by the four wind 
models (see Fig^), together with the interaction of the ex- 
ponential ejecta with a constant density ISM (triple-dot- 
dashed line). The ages of Tycho's SNR and SNR 1006 are 
marked above the horizontal axis. 
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Figure 5. Comparison of the ejecta density structure from 
different explosion models of SNIa with the exponential 
ejecta profile (triple-dot-dashed line), and a power-law 
profile with index 7 (long-dashed line), at time 5 x 10 8 s 
after the explosion. Each of the represented SNIa models 
is based on a different explosion mechanism: deflagration 
(dotted line), delayed detonation (solid line), pulsating 
delayed detonation (dashed line), and sub-Chandrasekhar 
(dash-dotted line). 



an exponential density profile for the ejecta. Instead, the 
simulated dynamics show a mixture of features from both 



these simpl e analytic ejecta profiles (see also Badenes & 



Bravo 2002|) 



Figure 4- The same as Fig.| but for the reverse shock. 



cannot be accurately described either by a power law or 



3. Discussion and conclusions 

We have shown that the CSM structure produced by dif- 
ferent presupernova wind models can be discriminated 
against an initially (at SN explosion) uniform ISM through 
the dynamical properties of the SNR. The expansion pa- 
rameter is most sensitive to the presupernova history at 
SNR ages larger than that of Tycho. We have also shown 
that different explosion models for SNIa can be discrimi- 
nated through the dynamical properties of the SNR, and 
through the different chemical abundances of freshly syn- 
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Figure 6. Expansion parameter of the forward shock as a 
function of time for the interaction of the analytical ejecta 
profiles with a constant density ISM. Due to the different 
ejected masses and kinetic energies of the models, the time 
scale has been normalized as in Dwarkadas & Chevalier 
(1998). The normalization timescale is 10.4 x 10 9 s. 
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Figure 8. Detail of the structures arising as a result of the 
interaction of a realistic explosion model (FDD has ben 
chosen for this example ) with a CSM modified by a time- 
dependent wind (wind model A here). The plot shows the 
density profile at the age of SNR 1006, in the region be- 
tween the ejecta reverse shock and the wind forward shock. 
The location of some secondary shocks has been marked 
with small arrows. 
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Figure 7. The same as Fig. \Q but for the numerical ejecta 
profiles (DEF, DDT, PDD, and SCH). Table | gives the 
corresponding normalization timescales. 

thesized elements (Fe, Ca and intermediate mass elements 
for X-ray lines, 44Ti for 7-ray lines). The expansion pa- 
rameter is most sensitive to the explosion mechanism for 
SNR ages shorter than that of Tycho. 

The interaction of realistic CSM, produced by non- 
constant presupernova mass loss, with the ejecta density 
profiles obtained from realistic SNIa hydrodynamic cal- 
culations results in SNR configurations that show a lot 
of structure. Particularly interesting is the presence of 
several secondary shock waves traveling outwards and in- 
wards (see Figure ||) . These secondary shocks could have 
important consequences for the overall SNR structure and 
properties (magnetic field, cosmic ray acceleration, hy- 
drodynamic instabilities, ionization structure, etc.). The 
denser areas in the hot, ionized intershock region will dis- 
play the highest X-ray luminosity (regions outlined by the 



dashed vertical lines). Thermal emission is expected from 
these regions, dominated by lines from solar abundance 
plasma for the shocked CSM (density spike close to the 
forward shock) or by lines from the elements synthesized 
in the SN explosion for the shocked ejecta (density spike 
between the contact discontinuity and the reverse shock) . 

Consistent modelling of the thermal X-ray emission 
from SNRs remains an unsolved problem. An accurate 
description of nonequilibrium ionization is necessary to 
obtain realistic tem perature and luminosity profiles (for a 
discussion, see e.g. |Borkowski, Lycrly fe Reynolds 2001 ). 
This is specially difficult for the heavy-element ejecta char- 
acteristic of SNIa. We are working on the coupling of hy- 
drodynamical models with atomic databases to be able to 
produce spectra models that can be used to fit observa- 
tions and thus take advantage of the high sensitivity and 
spectral resolution of XMM-Newton and Chandra. 
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